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Abstract 

Muons comprise an important contribution of the natural radiation dose in 
air (approx. 30 nSv/h of a total dose rate of 65-130 nSv/h), as well as in 
underground sites even when the flux and relative contribution are signif- 
icantly reduced. The flux of the muons observed in underground can be 
used as an estimator for the depth in mwe (meter water equivalent) of the 
underground site. The water equivalent depth is an important information 
to devise physics experiments feasible for a specific site. A mobile detector 
for performing measurements of the muon's flux was developed in IFIN-HH, 
Bucharest. Consisting of 2 scintillator plates (approx. 0.9 m 2 ) which measure 
in coincidence, the detector is installed on a van which facilitates measure- 
ments at different locations at surface or underground. The detector was 
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used to determine muon fluxes at different sites in Romania. In particular, 
data were taken and the values of meter water equivalents were assessed for 
several locations from the salt mine from Slanic Prahova, Romania. The 
measurements have been performed in 2 different galleries of the Slanic mine 
at different depths. In order to test the stability of the method, also measure- 
ments of the muon flux at surface at different elevations were performed. The 
results were compared with predictions of Monte-Carlo simulations using the 
CORSIKA and MUSIC codes. 
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1. Introduction 

The muons in the atmosphere originate from leptonic decays of pions and 
kaons generated by the high-energy collisions which cosmic rays experience 
with the atoms of the Earth's atmosphere. Muons are unstable particles 
which decay into electrons and positrons, respectively, with two correspond- 
ing neutrinos (electron (z/ e ) and muon (V^) neutrinos) with a life time (in there 
own reference system) of r M = 2.2/is. As leptons, muons are less affected by 
hadronic interactions and interact weakly with matter. They penetrate large 
thicknesses of matter before they are stopped and subsequently decay. Hence 
they are historically known as "penetrating component" of secondary cosmic 
rays, even detectable in deep underground sites. The cosmic ray muon flux, 
defined as the number of muons transversing a horizontal element of area per 



unit of time 



is of interest for various branches of science, in elementary 



particle physics as "heavy electron", as messenger of astrophysical processes, 
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in environmental and material research inducing natural radiation damages, 
and with a role for cosmogenic production of long living isotopes. The fo- 
cus of this paper is on studies of the atmospheric muon flux in the Slanic 
Prahova underground site. This site is actually under discussion as location 
of a large defeefor of the LACUNA p rOJ ee t flflfl. The pap, intends t o 
provide some basic information characterising the site, in particular on the 
underground depth of the salt mine. 

A different reason for measuring the flux of the atmospheric muons in 
underground arises from the practical necessity of information on the cosmic 
radiation background for different sites. This background does consists not 
only of muons which have survived the passage through the rock above, but 
also of contributions of natural radioactivity and of muon induced radiation 

ike neutrons, which can play a decisive role for low background experiments 

5]- 

For a simple and efficient procedure for measuring the muon flux at dif- 
ferent places a mobile device was set-up and operated since autum 2009, 
registering the muon flux at the surface and in the underground. Deter- 
minations of the water equivalent depth of any underground site could be 
done in a reasonable time scale. This feature is important in order to estab- 
lish very accurately the overburden thickness in water equivalent of matter 
(mwe). First measurements have been performed on the underground site 
of the Slanic-Prahova mine where IFIN-HH operates a low - radiation level 
laboratory []. 
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2. The LAGUNA project 



LAGUNA |2J (Design of a pan-European Infrastructure for Large Appa- 
ratus studying Grand Unification and Neutrino Astrophysics) is a research 
project , supported by the European Union to setup the infrastructure for a 
large underground laboratory with a first step to explore adequate locations 
for looking for extremely rare events like proton decay or for the experimental 
research of Dark Matter. 

Seven underground laboratories of Great Britain, France, Spain, Finland, 
Italy, Poland and Romania are involved. For LAGUNA, three detector types 
are considered based on different active detection media: MEMPHYS with 



water 
argon |9| 



, LENA, a liquid scintillator detector |8j and GLACIER using liquid 
. The site for LAGUNA experiments will be chosen along different 
criteria: the depth of the site i.e. the ability to absorb and shield against 
high energy muons, the available space and possibility to install a large vol- 
ume detector inside (larger than 100.000 m 3 ), and the natural radiation 
background. The site proposed in Romania is located in the salt mine Slanic 
Prahova, geographically situated at 45.23° N latitude and 25.94° E longitude. 
The elevation of Slanic site is 408 m above sea level at the entrance of the 
mine. The muon contribution to the natural radiation dose in air has been 
defined Q at ,1.3 ± 0.6 nSv/h at the Sam e e^tioo and .atitude as 
Slanic site. 

The salt ore from Slanic consists in a lens of 500 m thickness, few kilome- 
ters long and wide (see Fig. [1]). The salt is extracted from the Slanic mine 
continuously since ancient times and, due to this fact, many galleries (i.e. 
shaped caverns) are already excavated. 
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The largest one is the "tourist" mine "Unirea" (see Fig. [5]) characterised 

by: 

- depth: 208 m bellow ground level 

- temperature: 12.0 -13.0 °C 

- humidity: 65-70 % 

- excavated volume: 2.9 million m 3 

- floor area: 70000 m 2 

- average high: 52-57 m 

- distance between walls: 32-36 m 

- existing infrastructure: electricity, elevator, phone, Internet, GSM net- 
works. 

Besides the "tourist" mine UNIREA, muon flux measurements have been 
performed also in the "active" mine CANTACUZINO (see FigJT]) at two 
different levels with the physical depths of 188 and 210 m from the entrances. 

During the last 5 years a new laboratory [6[ for low background mea- 
surements was installed by IFIN-HH in the UNIREA salt mine of Slanic 
Prahova. From the Slanic site a huge volume of material has been already 
excavated, but the shallow depth could induce a problem. Following [10j for 
the GLACIER experiment the Slanic mine could be a feasible location as for 
this technique a depth of only 600 mwe is necessary. The main goal of this 
work was to determine the water equivalent depth (mwe) of the Unirea mine. 

3. Monte-Carlo simulation of the muon flux 

Monte-Carlo simulations were used to perform some preliminary explo- 
rations regarding the expected results of the experimental studies. Different 
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simulation codes have been used: 

- CORSIKA 6.735 12] (COsmic Ray Simulation for KAscade), a sophis- 
ticated Monte-Carlo code for simulations of the development of extensive air 
showers (EAS) in the atmosphere, has been used to estimate the muon flux at 
surface. As input for the simulations the primary cosmic ray spectrum ( for 
different primary masses ) in addition to a model for the hadronic interaction 
of the primaries with the Earths atmosphere is necessary (see below). 

- MUSIC [13I (MUon Simulation Code) is a simulation tool for 3 dimen- 
sional simulations of the muon propagation through rock. It takes into ac- 
count energy losses of muons by pair production, inelastic scattering, brems- 
strahlung and ionisation as well as the angular deflection by multiple scat- 
tering. The program uses the standard CERN library routines and random 
number generators. The validity of the code has been tested comparing the 
simulation's results with experimental data of the azimuthal distribution of 
single muon intensities in the underground Gran Sasso Laboratory for zenith 
angles up to 60° as measured by LVD experiment The results 13] show 
a good agreement between the data and the simulations (see Fig. [3]). 

- GEANT 3.21 [l5j, the detector simulation package from CERN has 
been used to simulate the interaction of the muons with the detector and for 
a proper calibration of the signal. Since the muon capture was not included in 
GEANT 3.21, a modified version of this code developed in context of charge 
ratio studies Q, H] has been used for this study. Since both hadronic 
interaction models implemented in GEANT 3.21 (GEISHA and FLUKA 
[19I ) give reasonable results for muon simulation, for this study, the default 
model: GEISHA has been used. 
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The simulation code CORSIKA has been originally designed for the four 
dimensional simulation of extensive air showers with primary energies around 
10 15 eV. The particle transport includes the particle ranges defined by the life 
time of the particle and its cross-section with air. The density profile of the 
atmosphere is handled as continuous function, thus not sampled in layers of 
constant density. Ionization losses, multiple scattering and the deflection in 
the local magnetic field are considered. The decay of particles is simulated 
in exact kinematics and the muon polarization is taken into account. In 
contrast to other air shower simulations tools, CORSIKA offers alternatively 
six different models for the description of the high energy hadronic interaction 
and three different models for the description of the low energy hadronic 
interaction. The threshold between the high and low energy models is set to 
E Lab = 80 GeV/n. 



The influence of the hadronic interaction mode 



s, implemented in COR- 



SIKA, on the muon flux has been investigated in [24j], and all of them are able 
to reproduce the experimental muon flux data, with some differences regard- 
ing the charge ratio of muons. From this reason, CORSIKA results given by 



DPMJET 



221 ] model where compared with the new QGSJ 



3T2 



251 ] ones and 



with experimental findings obtained by WILLI detector 32J. The WILLI 
calorimeter (see Fig. HJ-left), installed in IFIN-HH Bucharest If]], is oper- 
ated since several years for measuring charge ratio and flux of atmospheric 
muons at energies smaller than 1 GeV (see Fig. HJ-left), particularly exploring 
its directional dependence. The results obtained by WILLI are affected by 
statistical and systematic errors, the last ones due to the large geometric 
acceptance of the detector (see Fig. fright). The experimental data agree 
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with the simulation results (see Fig. E]) in the limit of the statistical and 
systematic errors. For higher energies, the simulations have been compared 
with the experimental data of BESS detector [23j. Fig. [5]- left compares the 
results obtained by semi-analytical approximation with those obtained with 
Monte-Carlo calculations and with WILLI and BESS data. The muon flux 

ae of Nash 



on surface can also be estimated by semi- analytical formu ae of Nash [20] 
and Gaisser [2l|. The validity of formula given by Gaisser 2l) is restricted 
to muon energies > 10 GeV due the fact that the muon's free decay is not 
considered. The formulation of Judge and Nash [20] tuned to experimental 
results, may be useful for the consideration of low energy muon flux. 
In the present work, the DPMJET 



22] model was used to describe the 



hadronic interaction at high energy and, the UrQMD [26] model for low 
energy processes. Based on the fact that the Monte-Carlo simulation code 
CORSIKA reproduces very well the experimental result (at low energies) its 
output is used as input for the underground simulation. There is some sensi- 
tivity of the underground muon flux predictions to the hadronic interaction 
models used in the simulations. 

The primary particle fluxes were introduced in the simulation's set-up 
using the results observed by the AMS prototype during a space shuttle 
mission [2^], using the spectrum and the ratio of proton and helium nuclei 
simulated in separate runs. Fig. [6] shows the fluxes of primary protons 



(center) and helium nuc 
borne experiments 28, 



ei (right) measured by different balloon and satellite 



29 



31] . that were launched in different region 



with various geomagnetic cut-off. In this study, the AMS data have been 
used, because it covers a large range of cutoff rigidities from the geomagnetic 
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equator (maximum) down to vertical cutoff rigidities less than 0.2 GV. The 
simulated muon flux generated by CORSIKA, have statistical errors (due to 
the limited number of considerate primary particles) and some systematic 
errors. The main sources of systematic errors are the primary spectrum 
and the hadronic interaction models. Other errors, due to particle decay or 
particle tracking are negligible compared to the other ones. 

FigfT] displays the muon fluxes at surface (CORSIKA simulation) and in 
underground, assuming, in a first approximation, a flat overburden over the 
observation level corresponding to an equivalent depth of 600 mwe (simulated 
by MUSIC using the muon flux from CORSIKA as input). The triangles 
represent only the muons from surface that managed to pass through the 
rock. The energy cut off for the surviving muons is estimated to be around 
150 GeV. By these simulation studies we estimate the expected muon rate 
at 600 mwe to about 10 muons/m 2 min. 



4. The apparatus 

The mobile detector was set-up in IFIN-HH and it consists of 2 detection 
modules. Each module is a scintillator plate (NE 114 type) of 0.9025 m 2 



and 3 cm thickness, see Fig. |H1 divided in 4 parts (0.475 x 0.475 m 



2^ 
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34j, readout by two photomultiplier tubes (EMI 9902 type) which receives 
the signal trough a wave length shifter (NE 174 A type). The modules 
are arranged one on the top of the other (at 8.5 cm distance), in order to 
identify the transversing muons as coincidence event. The signals from the 4 
photomultiplier tubes are two by two OR-ed (1 or 2) and (3 or 4) and then 
are putted in coincidence (see Fig. [9]) using a gate of 50 ns so no correction 
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due to the dead time of the detector is necessary. A counter module registers 
the coincidence events. 

The detector response is simulated by use of the GEANT 3.21 code. The 
interaction of the muon with the active detector material and the deposit of 
the energy in the scintillator plates are analysed. Fig. [TO] displays the energy 
deposit by muons (generated from CORSIKA output) in each detection plate. 
The signal threshold is set to 1/3 of the most probable energy deposit (6.3 
MeV), i.e. to 2.1 MeV [34|. The energy calibration has been performed by 
minimising the normalised spectra of the expected energy deposit and the 
anode signal (see Fig JTT]) . 

Preliminary tests had been performed on surface in order to check the 
counting rate for each individual scintillator and the rate of coincidences 
between two scintillators. The disposal of the scintillators is illustrated in 
FigfTJ] The results are displayed in TabsJT] and |2j The rate of singles is rela- 
tively high compared with the coincidences, due to the radiation background 
and electronic noise. 



Table 1: The counting rates for each scintillator plate 



Scintillator plate 


1 


2 


3 


4 


Counting rate 


515 


550 


541 


508 



The acceptance of the detector was also investigated considering the muon 
flux given by CORSIKA, taking in to account the fact that not all muons 
interact with the first layer, but manage to pass to the second one (see Fig. 
fl~3l) . Thus, based on GEANT simulations, which include the signal threshold, 
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Table 2: The rate of coincidences between planes 



Scintillator plates in coincidence 


1-2 


3-4 


1-3 


2-4 


1-4 


2-3 


Counting rate 


0.56 


0.69 


61.54 


60.82 


8.75 


9.10 



a correction factor of +9% has to be applied on the observed muon rate. The 
muon flux is given by: 

M = oti • a 2 ■ R (1) 

where: ol\ = 1.09 - the correction factor given by simulation, and a 2 — 1.11 - 
the correction factor due to the detector's surface and R - the counting rate. 

The detector is installed on a van (see Fig. [T4|) allowing to move quickly 
the system. The electric power for the entire system is supplied by a mobile 
electric generator of 1 kW power at 230 V AC or by a 12-230 V inverter of 
1 KW power which transform 12 V CC from the car's battery to 230 V AC. 

Using this mobile detector, measurements of the muon flux have been 
performed on different altitude levels and geographical locations at surface 
and on different mines of the salt ore of the Slanic site. 

5. Measurements and results 

The measurements of the muon flux in the underground have been per- 
formed at the Slanic site at 3 different locations: in Unirea salt mine (in 
IFIN-HH lab - see FigTS]) at 208 m below the entrance and in the active mine 
Cantacuzino, in 2 different levels, first one at 188 m and the second one at 
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210 m below the entrance. All runs were performed at approx. the same 
hour of the day (noon) in order to reduce the influence of the solar activity 
and atmospheric conditions. The acquisition time for each data set was one 
hour. In Cantacuzino mine, where an access tunnel is available, the measure- 
ments have been performed using the detector installed on the van. In Unirea 
mine, the detection modules were removed from the car and transported by 
an elevator to the observation level. The results of the three measurement 
campaigns are displayed in Tab. [3J 



Table 3: The muon flux data obtained in underground measurements. The errors are 
purely statistics. 



Location 


Depth 
(from surface) 


Muon flux 

(m _2 s _1 ) 


mwe depth 


Unirea mine 
Cantacuzino mine - Level 8 
Cantacuzino mine - Level 12 


- 208 m 

- 188 m 

- 210 m 


0.18 ± 0.01 
0.19 ± 0.02 
0.09 ± 0.01 


610 ± 11 
601 ± 21 
790 ± 29 



The variation of the muon flux as a function of the water equivalent depth 
is given by [351 ] : 



V(X) = A ■ (X /Xy ■ exp{-X/X ) 



(2) 



where: A = 0.03, X = 1470 m.w.e. and rj = 2.5. 

The difference in the muon flux measured at approximately identical phys- 
ical depths in Cantacuzino (-210 m) and in Unirea (-208 m) is associated to 
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the different thickness of salt rock above the detection place. Unirea mine 
is consisting of a huge cavity up to 57 m between the floor and the roof. 
In contrast, the Cantacuzino mine has a relative homogeneous rock massive 
above. The differences between the muon fluxes in Cantacuzino and Unirea, 
are caused also by the fact that the overburden of the mine is not flat, due 
to the irregular hill's shape above the entrance. The salt ore has a relatively 
regular shape (see Section [2]) and a homogeneous composition consisting of 
salt (NaCl > 98%) and different impurities (< 2%) [36|. Based on that, 
the m.w.e. has been estimated using only the muon flux values. The errors 
reported in Tabj3]are purely statistics. 



Table 4: Measured muon flux at different elevations and locations (the altitude was de- 
termined with a GPS system). The errors are statistics, including some systematic errors 
due to meteorological effects. 



Latitude (deg) 


Longitude (deg) 


Altitude (m a.s.l.) 


muon flux (m 2 s x ) 


45.29 


25.94 


655 ± 5 


146,74 ± 8,24 


45.28 


25.97 


588 ± 5 


145,30 ± 8,14 


45.24 


25.94 


408 ± 5 


143,24 ± 7,97 


44.32 


28.19 


70 ± 5 


128,05 ± 7,19 


44.40 


26.10 


64 ± 5 


122,28 ± 6,76 


44.36 


28.05 


7 ± 5 


119,07 ± 6,69 



Taking advantage of the mobility of the system, the measurements of 
the cosmic muon flux have been performed for many locations on surface 
at different geographical positions and different elevations, from sea level up 
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to 655 m. The results are in good agreement with measurements reported 
previously in ref. {3?]], that estimate a flux of 127 muon/m 2 s at 259 m a.s.l. 
(above sea level). The results are compiled in Tab JH and displayed in Fig. [15j 
During the campaigns at altitudes 70 m and 408 m the observation conditions 
were different (wind and low temperature) compared to the others, which led 
to different values and larger error bars. The meteorological conditions could 
affect the low energy muons based on the variation of the air's density. A 
systematic error of aprox. 3 % has been estimated based on WILLI measure- 



ments [38j performed over a full month (see Fig JTB"]) . The underground muon 
flux cannot be affected by weather conditions due to the fact that only high 
energetic muons could achieve the underground observation level. The effect 
of the air temperature from the mine, which is constant all over the year, on 
the muon flux is negligible. The influence of the meteorological conditions 
on the muon flux is presently under investigation using a portable weather 
station. 

6. Conclusions 



Suggested by the muon flux measurements reported for other sites |39| . 
the water equivalent depth of different places of the Slanic underground site 
were determined. The water equivalent depths of the Slanic mine are 610 
m.w.e. for Unirea mine, 601 for Cantacuzino mine (level 8) and 790 m.w.e. 
for Cantacuzino mine (level 12), respectively. The Slanic site is a feasible 
location for the GLACIER detector to be located in Unirea mine, with respect 
to the determined depth of 600 mwe (see Fig. [T71) . In addition, further 
promising locations for LAGUNA at the Slanic site are under consideration. 
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A new cavern, 100 m below the Cantacuzino mine (see Fig. [T]) could be 
excavated in a reasonable time scale. In this case a depth of about 1000 mwe 
would be at disposal for experiments. 

In near future, further measurements at different locations in Unirea mine 
will be performed, in order to get an improved overview on the variation of the 
Unirea mine's water equivalent depth. We expect that the muon flux varies 
for different locations of the mine due to the variation of the overburden at 
the Unirea mine. 

Such muon flux measurements could be also used for geological studies, 
e.g. to explore variations in the rock density above the observation level. The 
mobility of the detector implies a considerable practical flexibility of using 
the procedure of measuring muon flux differences for various aspects. 
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Figure 2: Map of Unirea mine, with the fiBq laboratory of IFIN-HH 
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Figure 3: Azimuthal distribution of single muon intensities in the underground Gran Sasso 
Laboratory for zenith angles up to 60 as measured by LVD (data points with error bars) 
and generated with MUSIC (dashed curve) [3] 
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Figure 4: The energetic (left) and angular acceptance (right) of the WILLI detector 
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Figure 5: Left: Compilation of the muon flux at surface: simulations, semi-analytical 
approaches and measurements; Right: Azimuthal variation of the muon flux measured by 
WILLI [if]] and compared with CORSIKA simulations (2 different hadronic interaction 
models). The error bars are statistical and systematic ones. 



21 




Figure 6: The primary cosmic rays spectra for protons and helium, used as input of 
the simulations (left). The fluxes of primary protons (center) and helium nuclei (right), 
measured in different balloon and satellite experiments. 
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Figure 7: The simulated muon flux at the surface and in the mine. 
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Figure 8: The detection module. Design of KASCADE 
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Figure 9: The electronic detection system. 
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Figure 10: The energy deposit in the scintillator plates (left - the top layer, right - 
bottom layer). 




24 



Figure 12: The disposal of the scintillators 




Figure 13: The schematically view of the muons interacting with the detection system, 
simulated with GEANT code (green lines are muons, red lines are secondary electrons) 
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14: Photo of the mobile detector mounted in a van. 
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Figure 15: Measured results of the muon flux variation with altitude in m a.s.l. The 
rectangle point represents the results from [37]. The error bars are the statistical and 
systematic errors. 
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Figure 16: Muon count rate as a function of day number of the year 2009 



38] divided into 



10 minutes intervals, for E > 0.4 GeV. Error bars represent statistical errors. 
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Figure 17: Muon flux as function of MWE depths for different underground sites. 
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